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Abstract 

Context. The white dwarf (WD) mass distribution of cataclysmic variables (CVs) has recently been found to dramatically 
disagree with the predictions of the standard CV formation model. The high mean WD mass among CVs is not imprinted 
in the currently observed sample of CV progenitors and cannot be attributed to selection effects. Two possibilities have 
been put forward to solve this issue: either the WD grows in mass during CV evolution, or in a significant fraction 
of cases, CV formation is preceded by a (short) phase of thermal time-scale mass transfer (TTMT) in which the WD 
gains a sufficient amount of mass. 

Aims. Here we investigate if and under which conditions a phase of TTMT before CV formation or mass growth in CVs 
can bring theoretical predictions and observations into agreement. 

Methods. We employed binary population synthesis models using the binary_c/nucsyn code to simulate the present 
intrinsic CV population. To that end we incorporated aspects specific to CV evolution such as an appropriate mass- 
radius relation of the donor star and a more detailed prescription for the critical mass ratio for dynamically unstable 
mass transfer. We have also implemented a previously suggested wind from the surface of the WD during TTMT and 
tested the idea of WD mass growth during the CV phase by arbitrarily changing the accretion efficiency. We compare 
the model predictions of the TTMT and the mass growth model with the characteristics of CVs derived from observed 
samples. 

Results. We find that mass growth of the WDs in CVs fails to reproduce the observed WD mass distribution. In the case 
of TTMT, we are able to produce a large number of massive WDs if we assume significant mass loss from the surface 
of the WD during the TTMT phase. However, the model still produces too many CVs with helium WDs. Moreover, 
the donor stars are evolved in many of these post-TTMT CVs, which contradicts the observations. 

Conclusions. We conclude that in our current framework of CV evolution neither TTMT nor WD mass growth can fully 
explain either the observed WD mass or the period distribution in CVs. 
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1. Introduction 

The class of compact binary stars comprises a great diver¬ 
sity of stellar objects and phenomena in the galactic zoo. 
They are very important probes of our comprehension of 
stellar evolution in general and mass transfer in particular. 
Cataclysmic variables (CVs) are compact binaries consist¬ 
ing of a white dwarf (WD) and a low-mass main sequence 
(MS) star that transfers mass to the WD by Roche-lobe 
overflow. CVs have been investigated for several decades, 
but their formation and evolution is still not fully under¬ 
stood. It is generally accepted that CVs result from wide 
binaries evolving into a common envelope (CE) structure, 
from which the core of the giant, that is, the primary, 
remains as a WD and in which the separation decreases 
sig nifican tly by means of drag forces within the envelope 
(lPaczvnskilll976D . After the envelope is expelled, the orbit 
of the detached post-common-envelope binary (PCEB) is 
further reduced through the loss of orbital angular momen¬ 
tum by gravitational radiation (GR) and magnetic breaking 
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(MB). When the orbit is sufficiently close, the accompany¬ 
ing MS star fills its Roche lobe, and if the resulting mass 
transfer is stable, a CV is born. 

According to this CV formation scenario, WDs in the 
newly formed CVs should have a mass distribution that 
is similar to the mass distribution of single WDs, if not 
shifted towards lower masses by an early expulsion of the 
envelope, which prematurely terminates the mass growth 
of the giant’s core. This naive expectation of on aver¬ 
age low WD masses ha s been confirme d by binary pop¬ 
ulation models of CVs, iPolitanol (ll996ll for instance pre¬ 
dicted a mean WD mass of 0.49 Mq for the primaries of 
CVs. However, measurements of WD masses in CV s have 
been in the range of [0.8-1.2 1 M 0 (e.g., I Warned ri97.lL 119761 : 
iRitterl 119761 : iRobinsonl I1976D , which is significantly higher 
than predicted. This discrepancy between the observed 
and expected mean WD mass in CVs has been success¬ 
fully i nterpreted as a selection effect bv iRitter &; Burkerf] 
(1986). Simply speaking, the idea is that the higher the 
WD mass, the more energy is released per accreted unit 
mass and the more extended is the accretion disk around 
the WD. Thus, CVs with massive WDs are (on average) 


1 





















T.P.G. Wijnen et al.: WD masses in CVs 


significan tly brighter and much easier to be discovered. 
However, IZorotovic et al.l (2011) recently showed that this 
previous explanation does no longer hold. They showed that 
the observed WD mass distribution of faint CVs (domi¬ 
nated by the emission from the WD instead of by that 
from the accretion disk) shou ld be biased towards low- 
mass WDs, while, as shown bv iLittlefair et ahl (l2008h and 
I Savoury et al.l (|2011h . the measured mean WD mass for 
these systems still remains at ~ 0.8 Mg. 

Thus the standard model of CV evolution might miss an 
important ingredient. IZorotovic et al.l (1201 lh suggested two 
possibilities. First, WDs in CVs may gain mass through ac¬ 
cretion of transferred matter over a nova cycle if less mass 
is expelled during the eruptions than is accreted between 
them. This contradicts standard theories for nova outburst s 

(]Prialnildl 1 986t iPrialnik fe Kovet3ll995t lYaron et al .1120051) . 

WD masses can only grow in CVs if no mixing of core 
matter w i th ac creted matter is assumed (IWilliamsI 2013; 
Starrficlcf 120151) . which is unrealistic. However, given the 
large discrepancy between the observed and predicted WD 
mass distributions, we need to investigate all possible sce¬ 
narios. Therefore, we here present binary population mod¬ 
els that include different mass accretion efficiencies in CVs. 


The second possi ble explanation put forward by 
IZorotovic et al.1 ( 2011 1 is that a large number of CVs 
descend from binaries with initially more massive sec¬ 
ondary stars. This implies a preliminary phase of ther¬ 
mal time-scale mass transfer (TTMT) in which the mass 
of the WD grows throu gh sta ble hydrogen burning on 
its surface (|Schenker et al.l [2002 .1. At that stage, the sys¬ 
tem might be observed as a s uper-soft X-ray source 
(iKahabka fe van den Heuvellll997il . A small sample of UV 
observations has shown that 10-15% of CVs accrete CNO 
processed material, indicating that the companion has been 
stripped of its external layers by a previous ph ase of TTMT 
(|Schenker fc Kindl2002t ICansicke et, al.ll2003ll . These com¬ 
panions therefore appear to be more evolved than a sin¬ 
gle MS star of the same mass. While binary population 
models of new-born CVs w ith evolved donor stars have 
been applied in the past (e .g. de Koo l 1992; Baraffc & Kolb 


[20001 [Podsiadlowski et all 20031 iKolb fe Willemsl 12005T ). a 


systematic study of the impact of TTMT on the WD mass 
distribution of CVs is missing. Here we fill this gap by using 
updated binary population models and investigate whether 
a large number of CVs descending from a phase of TTMT 
might explain the large masses of CV primaries. 


2. Code 

We simulated the evolution of a large number of binaries 
and selected those systems that evolve into a CV within the 
age of the Galaxy. Then we quantitatively investigated their 
characteristics and the evolutionary channel by which they 
have formed. We us ed the population nucleosynthes is code 
binary _c/nucsyn of llzzard et al.1 (120041 l200fil I2009H bas ed 
on the binary star evolution code of lHurlev et al.l (1 20021) . 

To simulate CVs, some specific aspects of their evolu¬ 
tion have to be taken into account. In particular, the imple¬ 
mentation and treatment of the stability of mass transfer, 
the mass-radius relation of the donor star, and the fate of 
the transferred matter have to be considered carefully. We 
address these issues in detail below. 



log M/M q 


Figure 1. Different regimes of stability against mass trans¬ 
fer for MS donor stars, depending on the mass ratio and the 
mass of the donor. The hatched region co rresponds to t he 
critical mass ratios for which, according to lPolitano (1996:), 
the mass transfer becomes dynamically unstable as a func¬ 
tion of the donor’s mass. The dotted line indicates the shape 
of this regions if one uses the condition that is originally in 
the code (see text). The solid curve represents the critical 
mass ratio for zero-age MS stars of solar composition above 
which mass transfer occurs on the thermal time-scale. This 
mass ratio is derived from the thermal mass-radius expo¬ 
nent by using the zer o-age MS mass-r adius relation imple¬ 
mented in the code ([Tout et al.l 1 19961) and assuming that 
mass transfer is conservative. 


2.1. Mass transfer 

The stability of mass transfer is determined by the change 
of the radius with respect to the Roche lobe. If the adiabatic 
response of the donor is unable to retain the star within its 
Roche lobe, mass transfer will occur in a dynamically un¬ 
stable way and will probably lead to a CE. However, if the 
star is able to restore hydrostatic equilibrium, mass transfer 
is determined by the thermal readjustments of the star. If 
the new thermal equilibrium radius exceeds the Roche lobe 
radius, mass transfer is driven by readjustments of the star 
on the thermal time-scale. Otherwise, the binary is stable 
against mass transfer and mass transfer is driven by angular 
momentum loss or nuclear evolution. 

The adiabatic mass-radius exponent for low-mass MS 
donors is very sensitive to the depth of their convective 
envelope. For low-mass MS stars (< O.7M 0 ) the envelope 
is deeply convective and the donor star is no longer able 
to restore hydrostatic equilibrium in response to mass loss, 
except for very low mass ratios. Therefore, the adiabatic 
mass-radius exponent and corresponding critical mass ra¬ 
tio decrease s teep l y around 0 .7 M e . FigurcQ] shows an an¬ 
alytic fit from lPolitanol (|1996f ) to the critical mass ratio for 
dynamically unstable mass transfer (hatched region): 

_ f | if M 2 < 0.4342 

qcl ~ (2.244(M 2 - 0.4342) 1 - 364 + § if 0.4342 < M 2 < 0.7 ’ 

( 1 ) 
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where M 2 represents the mass of the donor star in solar 
masses and the mass ratio is defined as q = M 2 . This fit 

is based on the adiabatic mass-radius e xponent as deter¬ 
mined by detailed model calculations of iHiellminel (Il989ll 
and is valid for conservative mass transfer. We used this 
more accurate prescription instead of the standard prescrip¬ 
tion in binary-c/nucsyn, that is, a constant q CI of 0.695 for 
M 2 < 0.7 Mg (which is also shown in Fig.|T] for compari¬ 
son). For masses above 0.7 Mq, mass transfer is assumed 
to be dynamic ally stable if q < 1.6 in our simulations 
(|de Mink et al.ll2007fl . This line is not shown in Fig.[l] 
Since binary_c/nucsyn does n ot follow the thermal 
mass-radius exponent of the donor dHurlev et al.ll2002l sec¬ 
tion 2.6.3) and mass transfer might not be conservative 
during TTMT, we used the mass-transfer rate to distin¬ 
guish between the CV phase and a phase of TTMT. We 
identify a binary as a system with TTMT if the primary 
is a WD, consisting of either helium (He), carbon-oxygen 
(C/O), or oxygen-neon (O/Ne), the donor is a MS star 
and the mass-transfer rate is higher th an the limit f or sta - 
ble hydrogen burning, as described in iMeng et al.l (2009). 
Furthermore, the mass of the WD has to increase by at least 
0.01 M 0 during the TTMT phase, otherwise we do not de¬ 
fine the emerging CV as a post TTMT system. Likewise, if 
the mass-transfer rate is below the limit for stable hydrogen 
burning, the primary is a WD, and the donor is a MS star, 
we identify the system as a CV. This definition of CVs may 
include some long orbital period systems with massive sec¬ 
ondaries (>0.7Mq) and lower mass WDs that just started 
thermally unstable mass transfer but did not yet reach the 
limit for stable hydrogen burning. We count these systems 
as CVs keeping in mind that they are apparently different 
from typically observed CVs. 


2.2. Mass-radius relation for CV donors 

Mass transfer can force the donor star out of thermal equi¬ 
librium when the time-scale of mass transfer is shorter 
than the thermal time-scale of the donor. In the case of 
low-mass donors, the mass transfer is driven by the loss 
of angular momentum due to MB and/or GR. We sub¬ 
tracted the angular momentum loss due to MB directly 
from the orbit, assuming th e orbit a nd spin are coupled. 
We used the prescription of iHurlev et al.l ([2002) for MB. 
MB is assumed t o be active until the d onor star becomes 
fully convective dRanoanort et, al.l 119831 . For single zero- 
age MS (ZAMS) stars and in detached binaries, this occurs 
when the mass is ~ 0.35 M 0 and MB is probably reduce d 
significantly in an abrupt manner (ISchreiber et al.l 1201011 . 
This disrupted MB scenario can explain the observed gap 
between ~ 2 and 3 hours in the orb ital period distribution 
of CVs (e.g. ISoruit fc Ritterl ll983l l. In semi-detached bi¬ 
naries that contain a secondary with a radiative core, the 
strong angular momentum loss and the resulting high mass 
transfer rates due to MB drive the MS donor out of thermal 
equilibrium, which causes the radius of the donor star to 
exceed its thermal equilibrium radius. As a result of this 
bloating and readjustments of the donor on relatively long 
thermal time-scales, the stellar structure of the donor cor¬ 
responds to the stellar structure o f a more mas sive single 
star on the MS (|Howell et al.ll200il lKniggell2006ll . In other 
words, the mass of the donor is lower than that of a MS 
star with the same radius and stellar structure. Therefore 


the mass-transferring donor in a CV becomes fully convec¬ 
tive at a lower mass than its MS counterparts in detached 
binaries or single stars. As a result, the dynamo mecha¬ 
nism responsible for MB rema ins activ e for dono r masses 
> [0.2 — 0.26] Mq (iMcDer mott fe Taamlll989tlHowell et al.l 

1200 ll ; iPatterson et al.ll2005l) . When MB becomes inefficient 
at an orbital period of ~ 3 hr, the donor has time to relax 
its radius to its equilibrium value, which is lower than the 
value of the Roche-lobe radius. The binary, having become 
a detached system, evolves towards shorter orbital periods 
driven by GR alone. At an orbital period of ~ 2 hr, mass 
transfer starts again, but at a much lower rate. In other 
words, the standard scenario for CV evolution explains the 
deficit of CVs in the period gap by predicting that CVs 
pass the gap as detached systems. 

Because of this standard theory of CV evolution, the 
mass-radius relation of the donor star is of crucial impor¬ 
tance for properly simulating the standard model of CV 
evolution, as the disrupted MB scenario only works if secon¬ 
daries above the gap have a larger radius than their MS ra¬ 
dius. Not assuming an increased radius above the gap would 
not allow simulating the orbital period gap seen in the ob¬ 
served distribution of CVs. Consequently, we would not be 
able to separate systems below and above the gap, which, 
as we show below, might be imperative for understanding 
the WD mass distribution in CVs. To account for the larger 
radius when the donor is out of thermal equilibrium, we im¬ 
plemented the mass-ra dius relation for low-mass MS donors 
in CVs, as deduced bv iKnigge et al.1 ( 2011 ). To establish a 
smooth transition between the equilibrium radius (i?2,eq) 
given by binary-c/nucsyn and the increased radius for CV 
donor stars (i? 2 ,cv)) we define the factor by which the ra¬ 
dius is increased with respect to its thermal equilibrium 
value, as a function of the current mass of the donor star 
m, 


/ M 


R2,cv(m) 

^2,eq(m) 


( 2 ) 


and let the radius grow exponentially with ti me towards 
the fully inflated value (i.e. R 2 ,cv) given bv IKnigge et al.l 
( 2011 ). The equilibrium radius excess (/ exc ) and current 
radius (-R 2 ) as a function of time and mass are thus given 
by 


/exc (tj m) = f(m ) + (1 - f{m))e 


(3) 


R 2 (t,m) = /exc(i, m)R 2 ,eq{rn), 


(4) 


where t is the time since the donor filled its Roche lobe 
and r is the time scale f or a ngular momentum loss in CVs, 
typically 10 7 years dDavis et al.1 200811 . I n acc ordance with 
the value for M conv from IKnigge et all (|2011h . we assume 
MB is disrupted at 0.20 M 0 for CVs. When MB becomes 
inefficient, the donor has time to relax its radius to its equi¬ 
librium value. In this case, we use Eqs. © and ©, replace 
i? 2 ,cv with the radius the star had just before i t detached 
and R 2 eq with the radius as described in IKnigge et al.l 
(|20H for donors below the period gap. This radius de¬ 
crease is implemented in the code in a similar fashion as 
the increase (Eq. [3]), that is, decreasing the radius expo¬ 
nentially, where now 


/ exc (f,m) = 1 - (1 - /(m))e *. (5) 
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We assume that donor stars whose initial mass is lower 
than 0.35 M 0 do not experience efficient MB, analogous to 
single ZAMS stars. Therefore, we did not inflate the radius 
for these donors, but only used t he radius for CV do nors 
below the gap. As pointed out bv lKnigge et al.l (1201 111 , the 
power-law approximation for the mass-radius relation of 
MS donors breaks down for masses < 0.05 M 0 . We therefore 
only considered CVs with a secondary more massive than 
0.05 M 0 . 


2.3. CVs with evolved donor stars 


While the assumed mass radius relation of CVs is reason¬ 
able for CVs with MS donor stars, we have to take into 
account that CVs from binaries with more massive secon¬ 
daries, such as those that passed through a phase of TTMT, 
may have significantly evolved donor stars. 

As has been shown bv iPodsiadlowski et al.1 (l2003l h the 
orbital period evolution of CVs with donor stars that have 
a central hydrogen fraction below 0.4 significantly differs 
from those with less evolved secondary stars. In partic¬ 
ular, CVs with evolved donors do not become fully con¬ 
vective at an orbital period of ~ 3 hr as “norma l” CVs 
do, but at possibly much shorte r periods (|Baraffe fe Kolbl 
l200f)t IPodsiadlowski et aD l2f)03li . In fact, they may even 
bounce back and evolve to longer periods before becom¬ 
ing fully convective. As a consequence, CVs with evolved 
donors pass the period gap as accreting systems with a 
bloated donor star that is significantly hotter than a typi¬ 
cal CV dono r at the same period. Such systems have indeed 
been found (iTlrorstensen et al.1120021 iLittlefair et ~ahll2006 t 

ei^ 


iThorstensenl 20131 : iRebassa-Mansergas et all 201411 . Thus, 

if many of them are produced in our simulations, then the 
predicted orbital period distribution could be significantly 
affected. Incorporating evolutionary tracks for CVs with 
evolved secondaries in our simulations is beyond the scope 
of this paper. Instead, we simply keep track of the central 
hydrogen fraction of the CV donors, derived from their ini¬ 
tial mass and current age, and determine the fraction of 
systems in our simulation that conta in evolved donors as 
defined bv iPodsiadlowski et al.l (l2003ll . that is, with donor 
stars that have a central hydrogen fraction below 0.4. 

To evaluate whether a given predicted fraction of 
evolved systems would disagree with the observations, we 
derived an upper limit from the observed period distribu¬ 
tion. Assuming that all observed CVs in the gap are CVs 
with evolved donors (this is a strict upper limit as some 
CVs with un-evolved donors are born in the gap) and tak¬ 
ing into account that the number of CVs at the lower edge 
of the gap in the observed distribution increases roughly b y 
a factor of 3-5 (jGansicke et al.l [20091 : iKnigge et al. 201 ill , 
we estimate an upper limit for the fraction of CVs with 
evolved donors (that evolve through the gap as accreting 
systems) of ~ 20 — 30 per cent. If this fraction is exceeded, 
one would expect more CVs in the gap and a smaller in¬ 
crease of the number of CVs at the lower edge of the gap 
than is observed. 


2.4. Modelling the mass-transfer rate 

Binary_c/nucsyn is a parametrized stellar evolution code 
and does not calculate the various mass-radius exponents 
that are discussed in Sect. O Instead the calculation of 


the rate at which mass is transferred to the accretor, M tr , 
contains a numerical factor to ensure that the mass-transfer 
rate i s (numerically) steady (see Eq. (59) from iHurlev et al.l 
l 2002i l. However, this factor is too small to let the star fol¬ 
low its Roche lobe within a few per cent during stable mass 
transfer when the radius is inflated. In other words: when 
the system is a CV, the mass-transfer rate is too low to 
be self-regulating. Furthermore, this parametrization of the 
mass transfer rate does not correctly model mass trans¬ 
fer on the thermal time-scale of the donor. We therefore 
adapted this factor for both cases separately as described 
below. 


2.4.1. Stable CV mass transfer 

To prevent the star from overfilling its Roche lobe by more 
than a few per cent, we made the mass transfer prescription 
more sensitive to changes in th e do nor radius by multiplying 
Eq. (59) from [Hurley et al.l (l2002f) with a factor of 1000 for 
stable mass transfer when the binary is a CV. 

For a numerically smooth transition, we increased the 
factor from 1 to 1000 on the same time-scale that we in¬ 
flated the radius of the donor star. This allows the mass- 
transfer rate to be self-regulating during the CV phase 
and prevents numerical instabilities in calculating the mass- 
transfer rate. 


2.4.2. Thermal time-scale mass transfer 

In the case of TTMT, the original prescription of 
iHurlev et al.l (12002 1 underestimates the mass transfer rate. 
We therefor e included the same factor as derived by 
IClaevs et al.l (2014), who have tested the resulting calcu¬ 
lation of the mass transfer rate against a deta i led binary 
stellar evolution code (ST ARS, lEggletonlll97ll : iPols et, all 
119951 : IcTebbeek et al.ll2008h . This provides a better descrip¬ 
tion of TTMT, although the mass transfer rate may be a 
factor of three higher than the maximum from the detailed 
STARS code, in which case the duration of TTMT is cor¬ 
respondingly shorter (see their Sect. 2.2.4). 

2.5. Response of the primary to mass transfer 

In the previous sections we have described the consequences 
of mass transfer for the mass-losing star. An equally crucial 
and precarious question is how the primary star reacts to 
accretion, in particular, if the mass transfer is conservative. 
In what follows, we describe our approach for both TTMT 
and mass transfer driven by angular momentum loss. 


2.5.1. Thermal time-scale mass transfer 

When the mass-transfer rate is within the limits of stable 
hydrogen burning, all transferred hydrogen-rich matter is 
processed into helium and accreted onto the WD. If the 
mass-transfer rate exceeds this limit, hydrogen will be ac¬ 
creted faster than it can be processed into helium. Two 
scenarios for this situation have been proposed in the liter¬ 
ature. The first is that the redundant hydrogen-rich matter 
accumulates onto the surface of t he W D and forms a red- 
giant-like envelope dNomoto et al.lll979ll . Consequently, the 
system will probably evolve into a second CE instead of be¬ 
coming a CV. The other scenario is that the burning of hy- 
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drogen may cause a very strong wind (iHachisu et al.lfl996i ). 
This wind ejects part of the accreted matter and tends to 
stabilize the mass accretion onto the WD, thus prevent¬ 
ing the formation of a new giant-like envelope. A certain 
amount of matter will still be accreted onto the surface 
of the WD, at a rate M acc , depending on the mass accu¬ 
mulation efficiency of hydrogen burning, rj h, and the mass 
accumulation efficiency for helium-shell flashes, ?yne- This 
rate can be expressed as 

M acc = 77h r?He Aftr, (6) 

where M tr is the rate at which mass is transferred from the 
donor to the WD. The efficiency parameters 7 /h and ?7He 
depend on M tr (iMeng et al .1 1 2009! ) . We refer to this wind as 
the Hachisu wind. 

A second problem concerning the response of the WD to 
mass accretion occurs if the WD is made of helium. The ac¬ 
cretion of helium can occur either directly or through stable 
hydrogen burni ng on the surface of the WD. Models from 
IWooslev et al.l (I 1986i ) have shown that a detonation due to 
the accretion of helium at a rate of 2 x 10~ 8 Mgyr _1 can 
occur wh en the star reache s 0.66Mg, while previous mod¬ 
els bv lNomoto fc Sugimotol (Il977tj found a limit of 0.78 M 0 
for the accretion of hydrogen at the same rate. Based on 
these predictions, the maximum mass of a helium WD (He 
WD) was set to 0.7 M^ in binary c/nuc syn (|Hurlev et al.1 
2002). However, IShen fc Bildstenl (120091) showed that all 
of the pre-1991 results are calculated with an erroneously 
high value of the conductive opacity and the resulting 
ignition masses are a factor of 2 too low. Furthermore, 
ISaio fc Nomotol (I1998D found that the accretion of he¬ 
lium onto a 0.4 M 0 WD at a rate of 10~ 7 M Q yr _1 and 
10~ 6 Mgyr _1 , that is, the rate of TTMT, induces a stably 
burning helium star. This implies that a phase of TTMT 
can eventually turn a low-mass He WD into a high-mass 
carbon-oxygen WD (C/O WD). Since there is currently no 
clear stringent limit on the mass of He WDs, if there is any, 
we did not assume a mass-limit for accreting He WDs. 


2.5.2. Mass transfer driven by angular momentum loss 


When the mass-transfer rate is below the limit of stable 
hydrogen burning, the accreted hydrogen is compressed 
onto the surface of the WD and is subsequently ignited 
under highly degenerate conditions. This leads to unsta¬ 
ble hydrogen-shell burning and flashes, that is, to nova 
eruptions. As mass transfer can be assumed to be stable 
and conservative in between two (virtually) instantaneous 
nova eruptions, the critical mass ratio given by Eq.|T| can 
be used to distinguish CVs from systems that are unsta¬ 
ble against dynamical mass transfer, and might evolve into 
a second CE. However, it has been subject of debate for 
several decades how much mass is lost during outbursts. 
The long-standing paradigm has been that nova eruptions 
in CVs expel most of the mass that has been acc reted be¬ 
fore the nova outbursts and p r obably even more dPria lnik 
H)8^_ |Prialni^fc Kovetd IT 9951 iTownslev fc Bildstenl 12004 


Yaron et al 


Yaron et al 


200517 We used m a ^ _... CJ _ 

(120051) to construct an interpolation table with 


and 7Ti e j from Table 2 in 


efficiencies for mass ‘accretion’ during nova cycles, which we 
implemented in binary-c/nucsyn. The amount of mass loss, 
or in other words, the ‘accretion’ efficiency, depends on the 
mass of the WD, the mass-transfer rate, and the core tem¬ 


perature of the WD. We derived the core tempera ture of 
the WD from its mass and luminosity (lMestellll952 ). 

As outlined in the introduction, the immense discrep¬ 
ancy between observed and predicted WD mass distribu¬ 
tions in CVs motivated us to investigate the implications 
of WD mass growth during nova cycles in a separate model. 
In these models we arbitrarily increased the mass accretion 
efficiency during a nova cycle from 10 to 100%. 


3. Binary population synthesis 


We generated a three-dimensional grid with Mi, M 2 and 
the separation as free initial parameters. Our resolution was 
150 for each parameter, thus for each model we simulated 
~ 3 • 10 6 binary systems. We let the initial mass of the pri¬ 
mary, Mi, range from 1 to 9 Mg, to let the primary evolve 
into a WD within the Hubble time. T he i nitia l mas s func¬ 
tion (IMF) of the primary is given bv lKrouna et al . (1199 30. 
For M 2 we assumed an initial mass-rati o distribution that 
is flat in mass ratio q (ISana et al.ll2009l ). The initial mass 
of M 2 ranged from 0.08 to 3.5 Mg. Both Mi and M 2 were 
picked from their initial mass distribution with a logarith¬ 
mic spacing. The i nitial orbital separation a was assumed 
to be flat in log a dPonova et al.lll98H iKouwenhoven et abl 
120071) and ranged from 3 to 10 4 Rq to cover the whole space 
of binaries that will interact within the Hubble time. Each 
binary in the grid was given an individual formation proba¬ 
bility, depending on the assumed distribution of the initial 
parameters and taking into account the logarithmic spac¬ 
ing in mass. We furthermore assumed circular orbits and 
solar metallicity for all binaries. We set the CE efficiency 
parameter, a ce , equal to 0.25, in accordance wi th the range 
of values determined bv IZorotovic et ah (l2010i). For details 
of t he CE prescription w e used, see iHurlev et al.1 (I2002D 
and IZorotovic et al.l (l2010h . The binaries were formed with 
a constant star formation rate by randomly assigning them 
a lifetim e betwee n 0 an d 13.5 Gyr, our assumed age of the 
Galaxy (iPascniini et al.ll2004l) . 


The motivation for this work is the disagreement be¬ 
tween observed_aiid predicted WD mass distributions of 
CVs. IZorotovic et al.l (l201lh suggested two possible solu¬ 
tions for the high WD masses in the observed sample, and 
we calculated three models to test both ideas against the 
standard scenario of CV evolution. In the first model, we as¬ 
sumed that a giant-like envelope is formed when the mass 
transfer rate exceeds the limit for stable hydrogen burn¬ 
ing, which should produce relatively few post-TTMT CVs. 
This model, which we refer to as the reference model, con¬ 
tains the most common assumptions of CV formation and 
evolution as discussed in the previous sections. In our sec¬ 
ond model, we added the stabilizing Hachisu wind dur¬ 
ing TTMT as this is known to be fundamental to produce 
large numbers of post-TTMT CVs, which potentially con¬ 
tain massive WDs. In a third model, we incorporated net 
mass growth during nova cycles into our reference model in¬ 
stead of the ‘accretion’ efficiencies derived from lYaron et al.l 
&■ 

The three main models can thus be summarized as fol¬ 
lows: 


1. Our reference model based on typical assumptions for 
CV evolution. 
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2. Our reference model, including a wind from the accret¬ 
ing WD that stabilizes mass transfer and prevents the 
formation of a giant-like envelope (Sect. 12. 5. ID . 

3. Our reference model, but assuming different efficiencies 
of mass accretion onto the WD. 

4. Results 

We first describe the results of our reference model for CV 
evolution. Subsequently, we address the effects of including 
the Hachisu wind, and then the mass growth. 

4.1. Reference model 

Figure^] shows the binaries that are identified as a CV at 
the present epoch (according to our definition of a CV, 
see Sect. ED- The distribution of both the orbital period 
(bottom) and WD mass (left) are shown, as well as their 
combined probability distribution in the two-dimensional 
plane. The vertical dashed lines mark the upper (3.18 hr) 
and lower edge (2.15 hr) of the observed period gap (|Kniggel 
120061 ). The hatched region indicates the parameter space in 
which mass transfer will become dynamically unstable, as¬ 
suming q CI as given as in Eq. ©• The dotted line marks the 
shape of the dynamically unstable region if one assumes a 
constant q CI of 0.695 for M 2 < 0.7 Mg. These regions are 
depicted in the same way as their corresponding regions 
in Fig© and were calcula ted by eq uating the radius of the 
donor star, as given bv lKnigge et all (1201 lh . with its Roclie- 
lobe radius and solving for the orbit with Keplers third law 
using the corresponding critical mass ratio. The right bor¬ 
der of this region corresponds to a donor mass of 0.7 M 0 , 
the curved left border corresponds to the critical mass ratio. 
Some CVs are located in the region for dynamically unsta¬ 
ble mass transfer. To explain this, we have also shown the 
right border for this region if one assu mes the mass-ra dius 
relation that is originally in the code (iTout et alJll996h . in¬ 
dicated by the dashed-dotted line. Since virtually all of the 
CVs in the hatched region are located on the stable side 
of the dashed-dotted border, it illustrates that these CVs 
have formed recently an d their radius is st ill inflating to¬ 
wards the value given bv lKnigge et al.l (1201111 . The mass of 
the CV donors in this region is higher than 0.7M o and the 
system is currently still stable against mass transfer. These 
CVs, however, will experience dynamically unstable mass 
transfer when the donor star becomes deeply convective at 
a donor mass of ~ O.7M 0 . 

The orbital period distribution shows a clear spike at the 
period minimum, corresponding to very low-mass secon¬ 
daries (M sec < 0.08 Mg). Most CVs in our standard model 
are currently below the period gap, but only 33.6 % of the 
current CV population was born in or below the gap. This 
is a well-known predi ction of binary population models of 
CVs (e.g. IKoTG fl993h : because angular momentum loss is 
assumed to be much lower below the gap than above the 
gap, the evolutionary time-scale is shorter above the gap 
than below by one to two orders of magnitude. 

The absence of CVs with low-mass WDs in the ~ 
3 — 4 hr orbital period range, that is, in the hatche d region 
in Fig© has been predicted by previous models (Ide Kooll 
1 9921 1 and can be explained as follows: mass transfer be¬ 
comes dynamically unstable for CVs whose mass ratio ex¬ 
ceeds the critical value given by Eq. ©. This will oc¬ 
cur within the region marked by the solid line in Fig. [2] 


CVs are thus not able to evolve towards shorter periods 
through this region without experiencing dynamically un¬ 
stable mass transfer and merging. This means CVs with a 
WD < 0.7M 0 and < —0.7 have to be born there 

or must have descended from higher WD masses. 

The WD mass distribution shows three peaks. The first 
peak at 0.4 M 0 consists of He WDs that evolved into a CE 
when the primary was on the first giant branch (FGB). 
The second peak at 0.55 M 0 corresponds to binaries for 
which the mass growth of their core was terminated when 
they were on the asymptotic giant branch (AGB). These 
are C/O WDs with a lower mass than expected for single 
WDs as a result of the earlier expulsion of the envelope. 
The steep decline after the peak for C/O WD s agrees with 
the o bserved mass distribution of single WDs dKenler et al.l 
l2007lb The third peak predicted at 0.75 M 0 is atypical for 
the mass distribution of single WDs. This peak contains a 
significant fraction (32.5%) of CVs that have experienced 
TTMT and thus originates from CVs that initially resided 
at longer orbital periods and that had more massive secon¬ 
daries. It is a direct consequence of using Eq. m for the 
critical mass ratio at which mass transfer becomes dynam¬ 
ically unstable instead of a constant value for q CI , as is as¬ 
sumed by default in binary-c/nucsyn. The region in which 
mass transfer becomes dynamically unstable for a constant 
q cr of 0.695 for M 2 < 0.7 M 0 extends towards high-mass 
WDs as indicated by the dotted line in Fig. [2] Thus, in this 
case, all CVs with a WD < 1M 0 and log(j^) > —0.6 
would merge eventually. The occurrence of this third peak 
has been previously predicted by CV population models, 
using t h e sa me prescription for the critical mass ratio (see 
lPolitanolll996l his Fig. 4). If this prescription had not been 
used, the predicted WD mass d istribution would be similar 
to the one derived bv Ide Kooil (Il992h . 

Our reference model predicts that 15% of the CVs 
have undergone TTMT before evolving into a CV, which 
is consistent with the results from UV observations from 
iGansicke et al.l (I2003H . The relatively small fraction of 
such post TTMT CVs is caused by the formation of a 
giant-like envelope if the mass-transfer rate of hydrogen- 
rich matter exceeds the limit for stable hydrogen burning. 
Consequently, He WD primaries evolve into FGB stars and 
C/O or O/Ne WDs primaries into AGB stars. Therefore, 
the binaries do not evolve into CVs, but eventually evolve 
into a second CE and merge. In other words, the limits for 
the mass-transfer rate, between which hydrogen burning 
on the surface of the WD is stable, define a range that is 
very narrow. Mass-transfer rates are thus typically either 
too high, causing the CVs to merge, or too low to allow the 
WD to gain a sufficient amount of mass. This limits the 
impact of post TTMT CVs on the WD mass distribution 
of CVs in our reference model. 

4.2. CV populations including the Hachisu wind 

In our second model the implementation of the Hachisu 
wind (see Sect. 1 2.5. Til prevents the WDs from forming a 
new giant-like envelope in the case of TTMT. This wind 
thus serves as a tool for investigating the consequences of 
extending the narrow range for stable hydrogen burning to 
higher mass-transfer rates, that is to say, extending the pa¬ 
rameter space for TTMT systems that may evolve into CVs. 
The resulting WD mass (left) and orbital period (bottom) 
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Figure 2. Centre: Two-dimensional histogram of the orbital period and WD mass distributions in CVs for our reference 
model. The colour intensity represents the sum of the formation probabilities of all CVs residing in that two-dim ension al 
bin. The vertical dashed lines mark the upper (3.18 hr) and lower edge (2.15 hr) of the observed period gap ((Kniggei 
120061 1. The hatched region indicates the parameter space in which mass transfer will become dynamically unstable, 
assuming q cr is given as in Eq. ©• The dotted line marks the shape of the dynamically unstable region if one assumes 
a constant q cr of 0.695 for M 2 < 0.7 Mg. These regions are depicted in the same way as their corresponding regions in 
Fig-El The dashed-dotted lin e mark s the region where mass transfer becomes dynamically unstable if one assumes the 
ZAMS mass-radius from I Tout et al.l (119961) . Bottom panel: orbital period distribution of the CVs in our reference model. 
Left panel: WD mass distribution of the CVs in our reference model. 


distributions are shown in Fig. [3] As in Fig.[2j the vertical 
dashed lines indicate the location of the observed orbital 
period gap, while the solid and dashed lines in the main 
grey-scale plot represent the region in which mass transfer 
becomes dynamically unstable according to the two pre¬ 
scriptions for q CT discussed previously. The dashed-dotted 
line marks the dynamically unstable region, assuming the 
mass-radius relation from lTout et al.l (119961 1. 

The WD mass distribution for M wc j < O.7M 0 is practi¬ 
cally the same as the distribution predicted by the reference 
model shown in Fig. [2] and so are their evolutionary paths. 
Assuming a wider range of mass-transfer rates that allow 
stable hydrogen burning without the formation of a new 
giant-like envelope, however, leads to the prediction of a 
much larger number of CVs with massive WDs (> O.7M 0 ), 
which in the reference model evolved into a giant-like star. 
This increases the total number of predicted CVs by 50% 
with respect to the reference model and now almost half of 
all CVs, 46 %, have undergone TTMT. The mean WD mass 
of post-TTMT CVs increased significantly during TTMT: 
from 0.56 Mq before the phase of TTMT, to 0.86 M 0 af¬ 


terwards. When only the CVs with a WD > 0.7 M 0 are 
considered, 74.8 % experienced TTMT. 

The impact of the description used for q CI on the pre¬ 
dicted distributions is even stronger than in the reference 
model. Most WDs with masses between 0.7 M 0 and 1M 0 
would not evolve into CVs if one assumes a constant q CI 
of 0.695, because they would experience dynamically un¬ 
stable mass transfer when they evolve towards the region 
bordered by the dotted line (see Fig. [3]). In that case, the 
only possibility to form WDs in CVs with a mass around 
0.8 M 0 below the gap are initial primary masses > 3M 0 , 
which have a very low formation probability compared to 
less massive stars. 

The predicted WD mass distribution is fairly broad 
above the gap and more strongly peaked below, while the 
mean WD mass is comparable for both distributions (0.76 
and 0.70M 0 ). This is caused by two effects. First, as most 
CVs with high-mass WDs below the gap descend from 
CVs above the gap (meaning that they are not born be¬ 
low the gap), the assumed mass loss during nova cycles as 
described in Sect. 12.5.^1 slightly erodes these massive WDs 
(0.8 — 1.2 M 0 ) on their way from long (~ 9.5 hours, i.e. 


7 















































T.P.G. Wijnen et al.: WD masses in CVs 



5e-06 


4e-06 


3e-06 "o 

o 
cr 
S3 
g; 

2e-06 ^ 


1e-06 


0 


Figure 3. Same as in Fig. [2] but for model 2, in which a wind from the WD that stabilizes mass transfer is included 
('Sect. ETSTl) . This wind allows more binaries that experience TTMT to evolve into CVs. The prescription for q CI is crucial 
for the subsequent CV evolution. 


^days) ~ —0.4) to short orbital periods. Second, CVs 
containing low-mass WDs (0.4—0.6 M 0 ) above the gap can¬ 
not avoid dynamically unstable mass transfer when evolv¬ 
ing towards shorter periods, that is, they enter the hatched 
region in Figs.[T]and[3] when the secondary develops a deep 
convective envelope (at a secondary mass of ~ O.7M 0 ) and 
merge during a second CE. Consequently, the distribution 
of WD masses above the gap shows a wider spread than 
the distribution of WD masses below the gap. 

4.3. Mass growth during nova cycles 

In the third model, we investigate the other possibility 
to produc e many CVs containin g high-mass WDs men¬ 
tioned by IZorotovic et al.l (12011( 1: we tested whether net 
mass growth of the WD in a CV evolving through a nova 
cycle can shift the low-mass WDs towards higher masses. 
We used the same assumptions as in our reference model, 
but instead of the standard assumption of slow erosion of 
the WD, we here tested four different accretion efficiencies. 
Figured] shows the WD mass distributions assuming that 
10, 20, 50, or 100% (from top to bottom) of the transferred 
mass is accreted, which means that it remains on the WD 
after a nova outburst. These WD mass distributions are 
normalized by the sum of the formation probabilities of all 
CVs in each model. Assuming an accretion efficiency of only 
10 % has little consequences for the predicted WD masses. 


The mean mass increases slightly, but the overall shape of 
the distribution with its three peaks remains the same (the 
peak at 0.55 Mg is smoothed out because of the bin size). 
This is simply because the low-mass WDs with relatively 
high-mass secondaries still evolve into the dynamically un¬ 
stable mass transfer region and those with low-mass sec¬ 
ondaries gain very little mass. The situation remains rela¬ 
tively similar for an accretion efficiency of 20 %, but changes 
significantly if 50% of the transferred mass is assumed to 
contribute to mass growth of the WD. In the latter case, 
many CVs with low-mass WDs and initially massive sec¬ 
ondaries can circumnavigate the dynamically unstable re¬ 
gion and evolve into CVs below the gap. In addition, many 
WDs gain a significant amount of mass. As a consequence, 
the peak located at 0.75 Mq in the distribution predicted 
by the reference model largely dominates the distribution 
and is moved towards higher masses (~ l.OM 0 ). This ef¬ 
fect becomes stronger if all the transferred mass is assumed 
to be accreted by the WD. In this case, many CVs evolve 
into SNIa candidates with WD masses >1.2 Mg. The lat¬ 
ter model is, of course, not a realistic one (because a sig¬ 
nificant fraction of mass is obviously expelled during nova 
eruptions), but it serves to illustrate the potential effect of 
mass growth on the WD mass distribution in CVs. 
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Figure 4. WD mass distributions for different fractions of 
mass growth during nova cycles, i.e. assuming that more 
mass is accreted between two nova outbursts than is ex¬ 
pelled during the eruption. The fractions are normalized to 
the total formation probability of all the CVs in each model. 
A higher accretion efficiency shifts the peak at 0.8 M 0 to 
higher masses. 


5. Discussion 

Table[l] shows the statistics for our different simulations. 
Our results can be summarized as follows: the reference 
model predicts the WD mass distribution to be dominated 
by low-mass WDs and a third peak at ~ 0.75 M 0 where 
most post-TTMT CVs reside. The fraction of CVs that ex¬ 
perienced a previous phase of TTMT increased by a factor 
of three when we incorporated the Hachisu wind, which 
makes the peak at ~ 0.75 M 0 more prominent. Finally, the 
models in which we allowed for mass growth during the CV 
phase predict relatively small changes for low accretion effi¬ 
ciencies, that is, a slight increase of the number of systems 
containing high-mass WDs. High fractions of mass growth 
predict CV populations dominated by systems containing 
very massive WDs (>1.OM 0 ). 

Below we first address to what extend our results de¬ 
pend on the assumptions and simplifications we made on 
our modelling approach, before we compare our model pre¬ 
dictions with the observations. 

5.1. Model uncertainties 

The model predictions are sensitive to several assumptions, 
for instance to the critical mass ratio for dynamically unsta¬ 
ble mass transfer, the Hachisu wind, the calculation of the 
TTMT rate, the mass limit for He WDs, the initial mass- 
ratio distribution, the strength of MB, and the efficiency of 
CE evolution. 

5.1.1. Critical mass ratio 

We demonstrated that the simulation of CV evol utio n is 
extrem ely sensitive to the analytic fit for q CI from lPolitanol 
(|l996f h This analytic fit allows a significant number of CVs 
with a massive WD to evolve towards shorter periods, while 


they would have been regarded as dynamically unstable if 
one would have used a constant q CI for all low-mass MS 
stars, that is, a rough cut-off at ~ 0.7 M 0 . The critical 
mass ratio, and the corresponding region in which mass 
transfer becomes dynamically unstable, functions as a kind 
of ‘road block’ that prevents CVs with low-mass WDs to 
evolve towards shorter orbital periods. Instead, dynamically 
unstable mass transfer leads most likely to a second CE and 
the stars probably merge. The value of q CI therefore has a 
huge influence on the formation of ~ 0.8 M 0 WDs. Even 
fine-tuning all other parameters towards producing many 
post TTMT CVs, such as the Hachisu wind, would predict 
few WDs around 0.8 M 0 if q CI were constant. It is therefore 
of crucial importance to use an accurate value of q CI instead 
of a crude approximation. 

Since the critical mass ratio depends on the adiabatic 
mass-radius exponent, it is mainly determined by the stellar 
structure of the donor star and whether mass transfer is 
conservative or not. To our knowledge, the prescription we 
use for q CI is the most accurate currently available, but it 
is still based on two assumptions that might be critical for 
the formation of CVs, that is, a MS structure of the donor 
and conservative mass transfer. 

First, the detailed structure of the secondary stars 
evolving from TTMT, especially the mass at which a con¬ 
vective envelope develops, is quite uncertain because these 
systems were initially more massive. These donors are prob¬ 
ably driven out of thermal equilibrium by mass transfer at 
a high rate and their structure might thus to some degree 
correspond to that of a more massive star. If this is the 
case, the mass at which a deep convective envelope devel¬ 
ops, that is, where the adiabatic mass-radius exponent in¬ 
creases steeply, could be lower than the canonical value of 
0.7M 0 . Thus, for a given WD mass, the mass ratio of a 
system evolving from long orbital periods to shorter orbital 
periods could be lower when facing the dynamically unsta¬ 
ble boundary. Therefore more CVs with less massive WDs 
(<O.75M 0 ) would evolve towards shorter periods without 
merging because the left-shifted boundary provides a higher 
q CI for a given secondary mass. This would concern CVs 
born at long orbital periods and post-TTMT systems. The 
peak of the WD mass distribution at 0.75 M 0 in the refer¬ 
ence model (and in the model including the Hachisu wind) 
would probably be shifted towards slightly lower masses. 
The impact of shifting the mass limit at which the donor 
stars develop a deep convective envelope towards lower 
masses on the resulting WD mass distribution would be 
even stronger when the masses of the WDs are assumed 
to grow during the CV phase, which would allow low-mass 
WDs at long orbital periods to circumnavigate dynamically 
unstable mass transfer if they accrete a sufficient amount 
of mass. 

Second, we assumed that stable mass transfer (driven 
by angular momentum loss) during the CV phase is con¬ 
servative. However, CVs experience nova eruptions during 
which both mass and angular momentum is lost. If the an¬ 
gular momentum taken away by the mass expelled during 
a nova outburst significantly exceeds the specific angular 
momentum of WD material, the critical mass ratio would 
be lower for a given value of the mass-radius exponent than 
it would be in the case of conservative mass transfer. This 
would cause stable mass transfer to occur only in systems 
with low mass ratios and may solve the problem of the 
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Table 1 . Statistics of our different models. Columns 2 to 6: the fraction of binaries in the grid that currently are a 
CV, the average WD mass of the CV population, the fraction of CVs that have a central hydrogen fraction < 0.4, the 
average WD mass at the beginning and at the end of the TTMT phase and the percentage of CVs that had a phase of 
TTMT. Assuming a ‘stabilizing’ wind from the WD increases the current CV population by 50 % and triples the number 
of TTMT descendants. On the other hand, increasing mass growth allows more CVs that emerge from a phase of TTMT 
with a (on average) lower WD mass to circumnavigate dynamically unstable mass transfer. 




(Mwd) 

(M 0 ) 


CVs With TTMT 

Model 

CVs In Grid 

(%) 

Evolved Donors 

(%) 

(MwD, prior) 

(M 0 ) 

(MwD, after) 

(M 0 ) 

% 

Reference 

0.40 

0.62 

23.38 

0.51 

0.82 

15.0 

Hachisu wind 

0.63 

0.71 

37.22 

0.56 

0.86 

46.0 

Mass growth 
10% 

0.40 

0.65 

19.96 

0.54 

0.93 

7.3 

20% 

0.47 

0.71 

24.82 

0.51 

0.96 

8.8 

50% 

0.70 

0.90 

35.41 

0.48 

1.09 

9.1 

100% 

0.63 

1.08 

29.96 

0.48 

1.34 

11.2 


large number of predicted CVs with low-mass WD. We will 
discuss this possibility in a follow-up paper. 

This reasoning and the sensitivity of the results to the 
critical mass ratio imply that we may underestimate the 
relative number of massive WDs in our models. The adia¬ 
batic mass-radius exponent as shown in Fig.[l] is currently 
being scrutinized (Nelemans, Webbink, private communi¬ 
cation). A more accurate value for mass-transferring stars 
based on detailed stellar models would drastically reduce 
the uncertainties in our model predictions. 

5.1.2. Hachisu wind 


estimate the mass transfer rate. We therefore also ran our 
simulations with the original prescription, which underesti¬ 
mates the mass transfer rate. The outcomes differ in which 
binaries actually become a CV after TTMT, but the overall 
distributions of both the WD mass and the orbital period 
look the same. We therefore argue that using a me thod that 
models TTMT more adequately (see e.g. IChen et all 1 2 D1 -ill 
would not (significantly) change the WD mass distribution. 
The period distribution, on the other hand, may look dif¬ 
ferent since we have a large number of evolved donor stars 
in our CV population. We address this question below (see 
also Sect. es. 


Since there is a consensus on the regime in which 
the mass-transfer rate enables stab le hydrogen burning 
(iNomoto et al.l 120071 : IWilliamsl 1 201 .'ili . the only option to 
increase the number of CVs that descend from TTMT is 
to assume non-conservative mass transfer. The most pop¬ 
ular mechanism for mass loss during TTMT is a strong 
wind that reduces the mass accretion rate and allows more 
systems to maintain stable hydrogen burning on the sur¬ 
face of the WD (lHachisu et alJ 199 0 i. In the framework 
of this wind model, the mass accumulation efficiency for 
hydrogen burning, t) h, and in particular the mass accu¬ 
mulation efficiency for helium shell flashes, during the 
wind determine how much the WD is growing and thus how 
fast the mass ratio changes as the mass of the secondary 
decreases. These efficiencies therefore define the evolution¬ 
ary course along which a system evolves towards lower sec¬ 
ondary masses, that is, to the left in Fig.Q] Higher efficien¬ 
cies would imply a steeper decline of the mass ratio and 
more CVs with a phase of TTMT could circumnavigate 
dynamically unstable mass transfer. In contrast, lower effi¬ 
ciencies imply that CVs with a phase of TTMT are more 
likely to run into dynamically unstable mass transfer. Thus, 
the relative number of post-TTMT CVs, and consequently 
also the resulting WD mass distribution, depend sensitively 
on the mass accumulation efficiencies during TTMT. 


5.1.3. Calculating the thermal time-scale mass transfer rate 

The calculation of the mass transfer rate in binary_c/nucsyn 
does not depend on the mass -radius exp one nt. In the case of 
TTMT the prescription from lClaevs et all (l2014lf may over¬ 


5.1.4. Accretion onto He WDs 

An important uncertainty that has not been tested by our 
model calculations is the maximum mass of He WDs that 
accrete hydrogen rich material (that might be fused into 
He on the surface of the WD). According to the literature, 
this value is highly uncertain, as outlined in Sect. l2.5.H If 
He WDs can be more massive than 0.7 Mg, they provide a 
serious contribution to the observed peak at 0.8 M 0 and si¬ 
multaneously reduce the relative number of WDs <0.5 M 0 . 
In the model including the Hachisu wind (see Fig. [3]) the 
contribution of He WDs that accreted a significant amount 
of mass during TTMT to the peak at ~ 0.8 M 0 reaches 
24.9%. Thus, if there exists a mass limit for accreting He 
WDs, the predicted peak in the WD mass distribution will 
be significantly less pronounced. 


5.1.5. Initial mass-ratio distribution 


One of the most important and critical ingredients of binary 
population models is the initial mass-ratio distribution. We 
assumed that all initial mass ratios are e qually probable. 
Although both early dPonova et al. 1 982li as well as more 
recent observational works (Raghavan e t al.ll2010h seem to 
support a mass-ratio distribution favouring equal-mass bi¬ 
naries, the general consensus is that the initial mass ratio 
distribution is flat, even for massiv e stars that appear to b e 
most biased towards equal masses (ISana et al.ll2009i 12012 ). 
An initial mass-ratio that is flat in q does not favour the 
formation of nearly equal-mass binaries, which is essential 
to produce large fractions of post-TTMT CVs. An initial 
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mass-ratio distribution that is proportional to q would thus 
significantly increase the fraction of CVs containing high- 
mass WDs in the resulting distributions. This would cause 
the peak at 0.75 Mq to dominate the peaks at 0.55 and 
0.4 Mq more strongly, as is observed. 


5.1.6. MB and CE efficiency 

We have assumed MB to be inefficient for fully convec¬ 
tive secondary stars. This, so named, disrupted MB sce¬ 
nario has been well establishe d during the past decade 
both in single stars (see e.g. [Reiners fc Mohantvl I2012L 
and references t h erein) and in close compact binary stars 
(ISchreiber et alJ l2010h . However, the strength of MB 
when the secondary star contains a radiative core is not 
well known. Curre nt prescr i ptions dif f er by severa l orders 
of m agnitude (ISchreiber fe Gansickel 120031: iKnigge et alJ 
l201ll l. Reducing the strength of MB, for instan ce by consid¬ 
ering the normalization factor implemented by iDavis et al.l 
(|2008f h would increase the evolutionary time-scale of de¬ 
tached PCEBs with secondary stars in the mass range 
~ 0.35 — 0.7 Mq and slows CV evolution towards shorter 
periods for systems with donor star masses in the range 
~ 0.2—0.7 Mq. This implies that CVs with high-mass WDs, 
that is, those systems that form the peak at ~ 0.8 Mq in the 
model including the Hachisu wind, would dominate the sys¬ 
tems that become CVs below the period gap less strongly 
at shorter periods. Furthermore, owing to the slower evolu¬ 
tion towards shorter orbital periods of CVs above the gap, 
our models would predict more systems in this period range 
for decreased MB. 

Equally uncertain as the strength of MB, but perhaps 
even more important for close compact binary evolution, is 
the efficiency of CE evolution. It seems that, at least for 
low-mass secondary stars (~ 0.1 — 0.5 Mq), a relatively low 
value of the CE efficien c y (~ 0.25 ) agrees best with the 
obser vations (IZor otovic et a.l.l l2010h and numerical simula¬ 
tions (iRicker fc Taamll2012il . However, CE evolution is ex¬ 
tremely uncertain mostly because of three reasons. First, 
we have currently no idea if this low value also holds for 
higher secondary masses, second, we do not know if per¬ 
haps a post-CE circumbinary disk efficiently extracts an¬ 
gular momentum from t he binary or bit of the emerging 
detached PCEB (Kashi & Sokeilf201lf ). and finally, it also 
remains an open discussion if perhaps the internal energy 
stored in the envelope such as the recombi nation energy 
contrib utes to the envelop e ejec tion process (IWebbinkll2008l : 
iRebassa-Mansergas et al.ll2012lf . 

However, the main effect that the CE efficiency would 
have on the predicted WD mass distributions is that higher 
values lead to the prediction of more CVs containing He 
WDs. This is simply because a significantly larger fraction 
of initially close binaries that start mass transfer when the 
primary star is on the first giant branch can sur vive CE evo¬ 
lution . This has recently been confirmed bv IZorotovic et akl 
•: 201 1:. who showed that increasing a ce from 0.25 to 1 re¬ 
sults in an increasing fraction of PCEBs with He WDs and 
an extension of the WD mass distribution towards lower 
masses that cannot be obtained with lower efficiencies. 
Therefore, using a higher value of a ce would increase the 
disagreement between our simulations and observations. 


5.2. Comparison with observations 

The model uncertainties presented in the previous section 
imply that performing binary population models of CVs 
has to be based on several assumptions, and changing these 
assumptions affects the predicted distributions. We have 
shown that mass growth during the CV phase can shift 
large fractions to higher masses but struggles to predict 
a broad distribution with the mean value around 0.8 Mq. 
We have also shown that the latter can naturally be ob¬ 
tained if a large fraction of CVs are assumed to descend 
from TTMT, but a large population of CVs with He WDs 
remains in all our models compared to the observations. 
Below we compare in detail the model predictions with the 
observations of CVs and their progenitors. Because the or¬ 
bital period is an observable that is much easier to measure 
than the stellar masses, we first inspect whether the orbital 
period distributions predicted by our models agree with the 
observations before we take a detailed look at the predicted 
WD mass distributions. 

5.2.1. Orbital period distribution 

In general, our simulations predict orbital period distribu- 
tions that resemb le those of previous models (iKolbl Il993l : 
iHowell et al.ll200TI) and agree well with the main features of 
the observed distribution. The predicted period minimum, 
the accumulation of systems at the period minimum, the 
orbital period gap, and the prediction of most CVs resid¬ 
ing below the gap all agree with the largest ho mogeneous 
samp le of observed CVs currently available (Gansicke et al.l 
1 2009) . This agreement is expected as the predicted pe¬ 
riod minimum and period gap are a direct consequence 
of the mass-r adius relation we used for CV donor stars 
(iKnigge et al.ll20lD . which has been designed to reproduce 
the period gap and the period minimum. 

However, our models disagree with the observed dis¬ 
tributions with respect to some details. For example, we 
predict a deficit of CVs between log{ g^) = —0.9 and 
—0.7 but a number of CVs at longer orbi tal periods , while 
the opposite seems to be observed (Gansicke et al.l 12009) . 
Furthermore, the accumulation of systems at the orbital pe¬ 
riod minimum seems to be more pronounced in our simula¬ 
tions than in the observed sample. These disagreements are 
caused by both model uncertainties as described in the pre¬ 
vious section and, probably more important, observational 
biases and selection effects. For instance, while the pre¬ 
dicted relative number of CVs above the gap sensitively de¬ 
pends on the strength of MB and the CE efficiency, the dis¬ 
agreement between the predicted number of these systems 
and the observations is most likely dominated by observa¬ 
tional biases. CVs at the upper edge of the period gap are 
particularly bright (high mass-transfer rates) while systems 
at the period minimum are particularly faint (low mass- 
transfer rates). Both the larger number of systems at the 
upper edge of the gap and the less pronounced peak at the 
orbital period minimum are therefore probably caused by 
the fact that bright systems are more likely to be observed. 
This becomes clear when comparing the orbital period dis¬ 
tribution of the complete sample of CVs with the most 
homogeneous sample identified so far (see iGansicke et all 
2009, their Fig. 2). Similar results are obt ained if observa¬ 
tional biases are taken into account. In fact. lPretorius et al.l 
( 2007 ) have shown that observational biases impose a large 
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Figure 5. Mass distribution of W Ds in PCEBs. From to y 
to bottom: observed sample used in lZorotovic et al.l (|2011l) . 
distribution predicted by our reference model, distribution 
derived from the reference model including a strong wind 
during TTMT, distribution predicted if it is assumed that 
20% of the transferred mass during a nova cycle remains on 
the WD. The fractions are normalized to the total forma¬ 
tion probability of all PCEBs in each model. The predicted 
PCEB distributions are identical (because mass growth 
during TTMT or during nova cycles does not affect the 
WD masses of the progenitors) and agree reasonably well 
with the observations. 


Figure 6. Mass distribution of WDs in CVs for the different 
models described in this work and the distribution derived 
from ob servations. From top t o bottom: observed sample 
used in IZorotovic et al.l (1201 lh with the black histogram 
representing a sub-sample for which the mass determina¬ 
tion is presumably more reliable, distribution predicted by 
our reference model, distribution derived from the reference 
model including a strong wind during TTMT, distribution 
predicted if it is assumed that 20% of the transferred mass 
during a nova cycle remains on the WD. The fractions are 
normalized to the total formation probability of all CVs in 
each model. 


discrepancy between the predicted and observed period dis¬ 
tribution of CVs. Their predicted period distribution looks 
very similar to the period distribution of our models (see 
their Fig. 11). 

5.2.2. WD mass distributions of PCEBs 

As shown bv IZorotovic et ah ( 2011 1. the high WD masses in 
CVs are not imprinted on the observed sample of PCEBs 
containing low-mass (M-dwarf) secondary stars. This led 
the authors to the conclusion that either the currently 
known pre-CV sample is not representative for CV progen¬ 
itors or the masses of the WDs in CVs grow during nova 
cycles. In this paper we tested both possibilities. Before 
we directly compare the predicted WD mass distributions 
for CVs, we first test whether the PCEB WD mass dis¬ 
tributions of our models agree with the observed one. The 
corresponding distributions are shown in Fig. [5] where we 
define a binary as a PCEB if it had a CE, currently does not 
have Roclie-lobe overflow, and consists of a WD (He, C/O, 
O/Ne) and a MS star. We used relatively large bins because 
the uncertainties in the masses derived from observations 
are significant. As expected, the WD mass distribution of 
PCEBs does not depend on our assumptions because both 
mass growth during nova cycles or during a phase of TTMT 
occur after the detached PCEB phase. Furthermore, the 
simulated WD mass distribution of PCEBs lo oks simil ar to 
the di stribution of the PCEB sample from IZorotovic et al.1 
( 2011 1. Both distributions have the same tendency to peak 
at 0.5 M© and show similar scatter around their mean value. 
However, our models fail to reproduce the extremely high- 


mass WDs (>0.8 Mq) that seem to be present in the ob¬ 
served distributions. This is probably due to uncertainties 
in the masses derived from observations. The WD masses 
shown in the top panel of Fig. [5] have been measured us¬ 
ing th e spectr al deco mposition/fitting method described in 
iRebassa-Mansergas et al.l (|2007H . We recently learned that 
more robust measurements of the WD mass using eclipse 
light curves indicate a much lower mass for those systems 
were the spectra indicate massive WDs (see iParsons et al.l 
1201.31 for one example). Therefore, the disagreement be¬ 
tween prediction and observation concerning the massive 
WDs in PCEBs is most likely caused by limitations of the 
spectral decomposition/fitting method. We also conclude 
from Fig. [5] that, according to all three types of models, 
either virtually all CVs with a WD more massive than 
~ 0.8 M© have grown in mass due to mass transfer in¬ 
stead of being born this massive, or the whole population 
of PCEBs is not representative of CV progenitors and only 
a small fraction of them (those with more massive WDs) 
may evolve into CVs. 

5.2.3. CV WD mass distribution 

Figureini shows the observed WD mass distribution and the 
distributions predicted by the different calculations that 
we performed, which are normalized by the sum of the for¬ 
mation probabilities of all CVs in each model. We used a 
larger bin size than in Figs. [2] and[5] to facilitate compar¬ 
ison with the observed distribution. Our reference model, 
which is mostly based on the current hypotheses of CV evo¬ 
lution, drastically disagrees with the observed distribution 
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(second and first panel from the top, respectively). While 
a small peak at 0.8 M 0 is predicted, the distribution still 
contains far too many low-mass WDs and therefore also 
peaks at O.4M 0 , which is not observed. The mean WD 
mass predicted by the reference model, that is, 0.62 Mq, 
clearly is too low when compared to the value derived from 
observations: (Mwn) = 0.83 ±0.02 (|Zorotovic et al. 11201 ill . 
As shown bv IZorotovic et all (1201 111 , the high WD masses 
derived from observations cannot be explained by obser¬ 
vational biases, and therefore it seems that an important 
ingredient of CV formation and/or evolution is missing in 
the standard model. 

On e of the possibilities proposed by IZorotovic et al.1 
(2011) is to assume that a large number of CVs descend 
from TTMT. Indeed, as shown in the third panel from the 
top in Fig.[6l the mean WD mass increases to (Mwd} = 
0.71 M 0 if a strong wind that extends the range of mass- 
transfer rates leading to stable hydrogen burning is assumed 
during TTMT. However, the distribution is still not in ac- 
cordance with the mea n mass derived from observations by 
IZorotovic et al.1 (|201lh i.e. (Mwd) = 0.83 ± 0.02. Most im¬ 
portantly, there is still a significant contribution (20.9%) 
of low-mass WDs, in particular He WDs, to the WD mas s 
distribution, which is not observed (IZorotovic et al.ll201ltl . 

The second option outlined by IZorotovic et al.1 (12011 1 
is mass growth during the CV phase, that is, assuming 
that during a nova cycle the WDs gain mass instead of 
slowly being eroded as predicted by the standard model. 
However, while this does allow increasing the mean WD 
mass, none of the distributions calculated for different frac¬ 
tions of mass growth is similar to the observed one. In 
Fig. [6] (bottom panel) we show the distribution assuming 
that 20 % of the transferred mass is not lost during a nova 
eruption. This shows that it is in principle possible to ob¬ 
tain a WD mass distribution of CVs that is dominated by 
massive WDs assuming mass growth during nova cycles, 
but the predicted distribution still disagrees with the ob¬ 
served one. The peak at O.8M 0 is less dominant and less 
broad than in the observed distribution, and the predicted 
distribution also shows a peak for He WDs at 0.4 M 0 , which 
is not observed. 

Based on Fig.[6l none of the two options that we tested 
seems to be the definitive answer for the WD mass distribu¬ 
tion in CVs, especially regarding the lack of systems with 
He WDs observed. However, we did not perform any fine- 
tuning but just used standard assumptions plus the wind 
proposed by Hachisu during TTMT or allowed for mass 
growth of the WD mass during the CV phase. As discussed 
in Sect. 15. 1.51 one could further reduce the He WD fraction 
for example by reducing the CE efficiency or the critical 
mass ratio for stable mass transfer in CVs. The latter op¬ 
tion will be discussed in a second paper. 


5.2.4. WD masses as a function of period 

IZorotovic et al.l (|2011f ) not only showed the general WD 
mass distribution of CVs derived from observations, but 
also the WD mass distributions of CVs above and below 
the gap. This division shows that the dispersion of WD 
masses above the gap is larger than of WD masses below 
the gap, which strongly peaks at 0.8 Mq, while the mean 
masses of both sub-samples are nearly identical (see Fig.|7| 
left). 


The best model assuming mass growth disagrees with 
both observed properties, the predicted mean mass in¬ 
creases towards shorter orbital periods and both distri¬ 
butions are equally peaked (see Fig.JTJ right). In contrast, 
our model incorporating the Hachisu wind that produces 
a large fraction of CVs descending from TTMT and which 
best reproduces the overall observed WD mass distribution, 
also shows a wide spread of WD masses above the gap and 
a concentration around O.8M 0 below the gap (see Fig. [3 
middle). This could indicate that a significant fraction of 
post-TTMT CVs might be the missing ingredient in our 
current understanding of CV evolution. However, as men¬ 
tioned in the last section, the latter model still produces far 
too many CVs containing He-core WDs. 


5.2.5. Fraction of CVs with evolved donors 


ISchenker et al.l (|2002f ) were the first to predict that the sec¬ 
ondary stars of CVs that went through a phase of TTMT 
can look different from the secondaries of CVs that did not. 
First, the secondaries of post-TTMT CVs descend from rel¬ 
atively massive stars and might therefore be significantly 
more evolved than genuine low-mass stars. This should 
cause them to have later spectral ty pes above the ga p and 
earlier spectral types below the gap dKolb fc Baraffelll999(1 . 
The first is apparently the case for two famous long orbital 
period CVs, AE Aqr and V1309 Ori. The latter has be¬ 
come the standard explanation for some CVs below t h e pe¬ 
riod g ap with very early sp e ctral types dThorstensen et al.l 
l2002t iLittlefair et al.1120061 : iThorstensea 20131) . Recently. 
such a system has been found well inside the period gap 
(iRebassa-Manser gas et a ll 2014 b in agreement with the 
predictions of lKolb k, Baraffel (I1999H that CVs with evolved 
donors should enter the period gap at shorter orbital peri¬ 
ods or not at all. 

As a seco nd possi bility for identifying post-TTMT CVs, 
ISchenker et alJ (l2002f) predicted unusually high N/C UV 
line ratios for low-mass secondaries if CNO processed ma¬ 
terial is brought to the surface by convection. Shortly af¬ 
ter this prediction had been made, an HST/STIS snapshot 
program of 31 CVs with strong emission lines revealed that 
indeed four systems ha d ext remely enhanced N/C UV line 
ratios dGansicke et al.l120031) . 

Thus, there is clear observational evidence for the exis¬ 
tence of post-TTMT CVs. However, the fraction of systems 
with extremely enhanced N/C UV line flux ratios seems to 
be significantly smaller than the fraction of post-TTMT 
CVs suggested by the model including the Hachisu wind. 
Thus the model represents a viable explanation if only a 
fraction of post-TTMT CVs shows enhanced N/C UV line 
ratios. 

A second problem for the Hachisu model is that the 
predicted fraction of CVs containing evolved donor stars 
increases with the increasing number of post-TTMT CVs 
(see Table |T]) . This number would further increase if for 
instance a lower CE efficiency was used to reduce the num¬ 
ber of CVs with He-core WDs (which currently is too large 
compared with the observations). In Sect. l2.3l we estimated 
that the fraction of CVs with evolved donors (that evolve 
through the gap as accreting systems) should not exceed 
~ 20 — 30 per cent because otherwise the period gap would 
be less pronounced than is observed. Our simulated CV 
population assuming the Hachisu wind predicts almost 40 
per cent of CVs with evolved donors, which exceeds this 
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Figure 7. WD mass distribution used bv IZorotovic et ail (1201 If) (left) of our model with the Hachisu wind (middle) and 
of our model including 20 per cent mass growth of the WD during a nova cycle (right) separated into the distribution of 
CVs in or below the gap (top panels ) and CVs above the gap (Bottom panels). The black histograms in the left figure 
represent a sub-sample of which the mass determination is presumably more reliable. The fractions in the middle and 
right panel are normalized to the total formation probability of all CVs in our models. The observations show a dispersion 
of WD masses above the gap, while the WD masses below the gap are more concentrated around 0.8 M 0 . Both features 
can also be seen in the model that includes the Hachisu wind (middle), but not in the 20 per cent mass growth model 
(right). 


limit. We therefore conclude that assuming a large number 
of CVs descending from TTMT can at least partly explain 
the high WD masses in CVs but unfortunately - gener¬ 
ates problems with the classical explanation of the orbital 
period gap. 


6. Conclusion 

We have performed binary population models of CVs 
with special emphasis on the predicted WD mass distri¬ 
butions with the aim to evaluate possible solutions for 
the problem of the high-mass WDs in CVs identified by 
IZorotovic et al.1 ( 201 1.1. We investigated the possibilities of 
WD mass growth during a preceding phase of thermal time- 
scale mass transfer and assuming net mass growth of WDs 
in CVs during nova cycles. In the latter case the predicted 
WD mass distributions drastically disagree with the ob¬ 
served one. In the first case, we find the best resemblance 
to the observed WD mass distribution by assuming a strong 
wind that increases the range of mass transfer rates, leading 
to stable hydrogen burning (as is frequently assu med in the 
single degenerate scenario for SN la, lHachisn et al.lIl996i l. 
This model also explains the observation that the disper¬ 
sion of WD masses is larger above than below the orbital 
period gap. However, this model still predicts a large pop¬ 
ulation of CVs with He WD primaries and evolved donor 
stars, which clearly contradicts the observations. We there¬ 
fore conclude that WD mass growth during nova cycles or 
during a preceding phase of thermal time-scale mass trans¬ 
fer cannot solve the discrepancy between the observed and 
predicted WD masses in CVs. 
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